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for of the entanglement loss reported previously. 

PACS numbers: 03.67.Mn, 03.65.Ud, 78.47.jp 



I. INTRODUCTION 

Entanglement is one of the key aspects distinguishing 
quantum from classical physics. Its fragileness, due to in- 
evitable coupling of a quantum system, such as qubits or 
photons, to a classical environment, however sets limits 
to its applicability to quantum information and quantum 
communication technologies. It is therefore very impor- 
tant to understand what is entanglement most vulnerable 
to and what processes can avert or undo entanglement 
loss. 

Qubits are the fundamental building blocks of quan- 
tum information science, and the last two decades 
marked great developments both theoretically as well as 
experimentally [1, 2]. After many successes in that field, 
the current research frontier is the qubit-qubit entangle- 
ment, which requires either a direct coupling between 
qubits or an indirect one through an auxiliary system, 
for example a resonator [3, 4]. As a result of that the 
qubits can potentially entangle to, or disentangle from, 
each other depending on the system design, its parame- 
ters, interaction time, or even the type of environments 

: [5-18]. 

In this article we focus on the process of entanglement 
revival in a system of two qubits driven by a strong exter- 
] nal, classical ac field and simultaneously coupled a quan- 
tum resonator which indirectly couples the qubits. We 
previously discussed this phenomenon in Ref. [15, 16] and 
demonstrated that entanglement does not only need to 
decay as the system evolves in time, but the system can 
also periodically regain some of its initial entanglement. 
Previously, it was shown for this system that the disen- 
tanglement between the qubits may be a consequence of 
the cavity dissipation [17, 18]. In this manuscript, we 
demonstrate that this is is not the only mechanism lead- 
ing to entanglement loss. Specifically, we show that the 
mere presence of a qubit-cavity coupling results in dis- 
entanglement in the subspace spanned by the qubits and 
that the further coupling of the cavity to the electromag- 
netic continuum leads to an overall tripartite entangle- 
ment decay. Therefore the best way to understand the 
qubit entanglement dynamics is by looking at the phe- 




FIG. 1: (Color online). Schematic representation of the setup. 



nomena from a larger, multipartite perspective. With 
this work we aim at providing insight into entanglement 
transfer back and forth within a multipartite system sub- 
ject to dissipation. 

This article is structured as follows. In Sec. II, we 
quantitatively introduce the system of strongly driven 
qubits. We derive the equations of motion and present 
their solutions. In Sec. Ill, we outline and discuss entan- 
glement measures applicable to the tripartite analysis. 
In Sec. IV we quantify entanglement between individual 
subsystems in a dissipationless regime. Subsequently, in 
Sec. V, we study the effects of an imperfect cavity on en- 
tanglement formation, revival and loss among different 
subsystems. We close the article with conclusions. 



II. THE MODEL AND ITS DYNAMICS 

We consider the system of two identical qubits exter- 
nally driven by a classical field of the amplitude A and 
the frequency uJc strongly coupled to a single mode res- 
onator [19] (see Figure 1). This system is described by 
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the following Hamiltonian, 



qubits [comment 1], 



H = Ho + + Hj , 

Hn 
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where represents the level spacing of each qubit, w 

is the frequency of the resonator cigcnmodc, and gj is 
the coupling strength of the Jaynes-Cummings type in- 
teraction between the j*^ qubit and the eigenmode of 
the resonator. Additionally, we use cr^ (a^,a) to de- 
note the qubit (resonator) raising and lowering operators. 
Throughout this article we chose units where = 1 . 

We work under the assumption that the qubits are 
driven very strongly (ensuring their greater resistance to 
decay 7), and are moderately coupled to the cavity mode, 
i.e. ^ I:^ w, Wc 5 7- Moreover, we consider the cav- 
ity dissipation rate k to be the dominant source of deco- 
herence in the system. A realisation of these conditions 
can be found for instance with superconducting qubits 
where a; w 5 GHz, g w 100 MHz, and 7 « 1 MHz [20]. 

The Hamiltonian (1) is time dependent. To sup- 
press the time dependence, one can apply a number 
of unitary transformations [19]. First we go to the 
frame oscillating with the driving field frequency lOc using 

U = eyiY> (^—iujct {i-'^a + Yl,j'^j /'^^ ^iid further to the 

interaction picture (IP) V = e-'^"''+"''^*He'^"''+"''^* . 
Upon ignoring the quickly rotating terms cx e*^*'^* (the 
strong driving regime) the effective IP Hamiltonian be- 
comes (see Ref. 15 for technical details) 



(2) 



where we set Q = uj,, and 5 — lo — ujc- What we however 
see is that in the strong driving regime the coupling to 
the resonator no longer mediates the interaction between 
qubits, but is simply reduced to the qubit state depen- 
dent bosonic displacement generator, i.e. the state of the 
multi-qubit system determines the state of the resonator, 
but the state of the resonator never affects the state of 
the qubits. 

We model the evolution of the system with the Lind- 
blad type master equation 



p = C{p) = -i[V,p] + KV{p) 



(3) 



where T) (p) = apd^ — ^{a^d, p} is a Markovian dissipation 

operator. The system is initially in a direct product state 
of a cavity field coherent state |a) and a Bell state of the 



* = ^(l+ + ) + |--» 



(4) 



where we used the diagonal basis |±) = ■^(|e)±|£f)), 

which are the cigenstates of the Pauli ctx matrix. The 
solution to the Lindblad equation (3) is given in terms 
of the density operator pija-.kifs, where the Latin indices 
i and k {j and I) stand for the |± ) state of the first 
(second) qubit and the Greek indices indicate the state 
of the cavity. The solutions to (3) are found in Refs. 
[15, 16]. The only non-zero entries of the density matrix 
Pija;kip with the initial condition given by (4) read 
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where ar,i are the real and imaginary parts of the initial 
state of the cavity and g = gi+g2 is the effective qubits- 
cavity coupling. 

The coherent state considered here has a continuous, 
time dependent amplitude. Such state is represented by 
a vector spanning the whole of the infinite Fock space, 
making the qubits-cavity system 2 x 2 x 00 dimensional 
rendering some of the entanglement measures inapplica- 
ble. These different coherent states however can be writ- 
ten in bases foimd in [21]. Using the fact that every co- 
herent state is a single-parameter state, we can recast the 
two coherent states \a±) in a two-dimensional form by 
means of orthogonalisation through the Gram-Schmidt 
process. 



t) = K) 
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X{t) = {a+\a-) 



(7) 
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such that ( t I -I- ) =0 and where the inverse transforma- 
tion reads 



\a+ ) 



t) 



|a_) = Vl-Ixni>+X|t). 



(8) 



As a result the system is now reduced to 2x2x2 = 8 
dimensions. In the last subspace the bases definitions are 
time dependent, but the resulting set of states {| t ), | J, )} 
is orthogonal at any point in time. In this form, we can 
easily use the established entanglement formalism. 

The full system time-dependent density operator p in 
the X X space spanned by the qubits diagonal 
basis and the orthogonalised coherent state basis reads 



(9) 
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which is a sparse 8x8 matrix where the only non zero 
elements are contained in the 2x2 blocks 
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III. ENTANGLEMENT MEASURES 

It is not always easy to establish non-separability (en- 
tanglement) of a system based on the form of a density 
operator. The factors that play a role here are, among 
others, bases choice and dimensionality of the system. 
Fortunately, the last two decades brought developments 
in the field of entanglement measures [22-25]. There it 
was shown that stepping beyond 2x2 physical systems, 
where one oftentimes uses concurrence;, one has to ac- 
count for a greater number of correlations between indi- 
vidual players in a multipartite physical system [26] and 
can choose from entanglement witnesses, negativity, or 
the three-tangle. We devote this section to briefly review 
some of the most important aspects which we will find 
useful in out subsequent analysis. 

One of the first entanglement measures to be intro- 
duced and since then widely used for 2 x 2 dimensional 
systems is concurrence [22] . Its mathematical form given 

by 



C = max(o,V%"- V^- V^- V^) , (10) 

where are the eigenvalues of 

R = {ay ® <Jy) p* {oy ® (Ty) p, wlth Ai being the largest 



of them, (Ty being the Pauli y-matrix. The value of C 
ranges from zero (no entanglement) to one (maximum 
entanglement). This measure however no longer suffices 
when dealing with systems involving more than two 
two-dimensional subsystems. 

In order to study entanglement in the tripartite sys- 
tem, we use Horodeckis' separability criterion [25] and 
stemming from it negativity [24] to quantify tripartite 
entanglement. Using the partial transposition (in the 
second subspace) defined by 



p = '^oiijki\i){j\ 



ijkl 
ijkl 

this criterion states that the density operator of an en- 
tangled state upon transposition in one of the subspaces 
will have at least one negative eigenvalue. Negativity is 
then the sum of absolute values of negative eigenvalues 
ofpf^. 

Thus when studying a tripartite system composed of 
three subsystems A, B and C (in this case A and B are 
the qubits and C is the cavity, but the labeling is com- 
pletely arbitrary), we can find the degree of entanglement 
between the combined bipartite subsystem AB and sub- 
system C, by partial transposing the density operator 
PABC of the system in the basis states that span the sub- 
system C, and later adding up all of the absolute values 
of the negative eigenvalues of p^^c- ^ result we ob- 
tain negativity Neg(A_B|C), which when equal to zero 
corresponds to no (or bound i.e. a state with zero nega- 
tivity that is not separable) entanglement and when equal 
to ^ indicates maximum bipartite entanglement. To get 
the full picture of tripartite entanglement in this system 
we need to also calculate Neg{AC\B) and Neg{BC\A), 
where the partial transposition is made in the subsystem 
B and A basis respectively. 

Since the dimension of this system is larger than six, 
i.e. the limit imposed by the Horodeskis' separability 
criterion, we could encounter bound entanglement . We 
could avoid this subtlety by creating a map (g) — >■ 

which maps the entangled state | + + ) + | ) onto a 

superposition | -|- ) -|- | — ) reducing the dimensionality of 
the system, by removal of permanently empty rows and 
columns of the density operator. This will however prove 
to be unnecessary, as we will see later, that the only time 
when negativity is strictly zero is at the expected peri- 
odically distributed points in time 6t = 2Tm for integer 
n, when the qubits are completely disentangled from the 
cavity (see Figure 2 and Eq. 5); something that can be 
easily seen without invoking any entanglement measures 
formalism. Thus the excessive dimensionality of our sys- 
tem posses no problems with regards to using negativity 
as an entanglement measure. 

One drawback of the negativity is that it only provides 
information about entanglement of two parts of the sys- 
tem under partitioning of our choice and does not tell us 
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anything about the total entanglement present. Adapt- 
ing the approach of [27] we can use the sum of the bipar- 
tite entanglements 

TE = Neg (Q1CIQ2) + Neg iQ2C\Qi) + Neg iQiQ2\C) . 

(11) 

where we replace the arithmetic mean by a direct sum so 
that we can use TE = 1 as an easier reference point for 
how much entanglement was there initially in the system. 
As a result of this definition we get that < TE < 
|, where the lower bound indicates no and the upper 
bound indicates maximal tripartite entanglement, that 
of for example the GHZ^ state 

\GHZs) = ^{\000) + \lll)) , 

where all negativities Neg (Qi(52|Q3) = 
Neg(QiQ3lQ2) = Neg(Q2Q3lQi) = i 

The GHZ3 state shows a feature that will be important 
to our further discussion, namely tripartite entanglement 
sharing. In this state (as opposed to the VF-state) the 
individual subsystems share bipartite entanglement but 
partial tracing over one of the subsystem (loosing a qubit) 
results in a statistically mixed state (the W state results 
in a Bell state). 

This result is known as the monogamy of entanglement 
which states that a subsystem A maximally entangled to 
second subsystem B cannot simultaneously be entangled 
with another subsystem C. This has been first formu- 
lated by the CofFman, Kundu and Wootters [23] in terms 
of the inequality 

^{A\BC) ^ ^(AB) +^(AC) ' (12) 

where are the tangles (concurrences squared). Here 
C(A|BC) is found by reducing the dimensionality of the 
density operator p down to the subspacc spanned by the 
two eigenvectors of p' with non-zero eigenvalues of the 
BC subspace and C(^b) and Cf^AC) concurrences of 
bipartite subsystem obtained by partial tracing the total 
tripartite system over the subsystem C and B respec- 
tively. 

The inequality (12) can be used to define a three-tangle 
given by the inequality mismatch 

TABC = C'^A\BC) ~ ^{AB) ~ ^{AC) ■ 

This new quantity tells us how much of residual tripartite 
entanglement is there when all of the bipartite contribu- 
tions are taken away. It is easy to see from the definition 
of concurrence that the three-tangle ranges from zero (no 
shared entanglement) to 1 (completely inseparable state 
of GHZ3 type). 

From the solutions to the equation (3) we see that the 
qubits entangle with the cavity, which for S, k ^ and 
t ^ 00 leads to a perfectly entangled GHZ-like state, as 



the coherent state amplitudes undergo a shift in opposite 

directions, 

\^,a){0) = \^)(g>\ao), 

|*,a)(0 = ■^{\ + + )\ao-2igt) + \--)\ao+2igt)) 

^l={\++,A+) + \--,A.)) . 

This state is different from the GHZ-state since 
( j4+|^_ ) 7^ 0. Upon taking the partial trace of p{t) 
over the cavity, we observe that the diagonal entries of 
the density operator Trc(|^)(^'|) are unchanged, but 

the entries | + + ) ( | acquire time dependence e~^^ * 

which mimics the dephasing of the two qubit state. This 
is because as t grows the state of the system more and 
more closely resembles the GHZ-state, and taking the 
trace leaves the state in a completely mixed state to a 
larger extent. It is a continuous in time analogue of the 
GHZ state formation from the original Bell state. 

The effect of entanglement revival in this system is 
brought about by the presence detuning between the cav- 
ity and the resonantly driven qubits. Since we would be 
interested in periodic revival of entanglement for the re- 
mainder of out analysis we have to keep S ^ 0. In what 
follows we will mainly focus on the dissipationless case 
as it provides a very good insight into qualitative as well 
as quantitative aspects of qubits-cavity entanglement dy- 
namics. Later we will study the effect a combination of 
dissipation and detuning on the inter-qubit as well as the 
qubits-cavity entanglement. 

IV. DISSIPATIONLESS CAVITIES 

In the closed system (in the k — > limit ), the non- 
resonant interaction between the qubits and the cavity 
will result in formation of a coherent state with an am- 
plitude oscillating in time with frequency S. Under these 
conditions the complete state of the system is still repre- 
sented by Eq. 5 where the previously defined expressions 
are replaced by 

hi = 0, (13) 

_ gg(cos<5t-l)(4|-2ifc+2ia(sin(5t)cos(5t-l)/5 

where a and b are the real and imaginary parts of the 
initial coherent state amplitude, and the value of /12 (t) 
will have no effect on the result. Upon partial transposing 
expression (9) with respect to the cavity subspace we get 




where there is only one negative eigenvalue, and the neg- 
ativity takes the form 

Neg(QiQ2|C) = -(^l-e ^ j . (14) 
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FIG. 2: (Color Online) Negativities: qubits vs cavity 
Neg (Qi(52|C) for a qubits initiated in the \9 ) state (4). 
Plots made for gi = 92 ^ lOOMHz and 5 = {0.1,0.3,0.5, 1} 
GHz in blue, green, black, red (increasing dashing frequency) 
respectively. The peak of the maximum is always where 
5t takes values of odd integers of vr and their value is 



16(81 + 82)-' 



Taking the partial transposes in the individual qubit 
spaces, we get a symmetric result 



PQi 



Tr(p)^ = |<i>)($| f^PQ.c 
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_ 1 I 

~ 2 
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where Tr(p)^ denotes a partial trace over the cavity 
states and e{t) = (a+ja- ). This has to do only with 
the fact that the Tr )^ entries carry time depen- 
dence, while the populations i.e. the Tr (p+_|_;++)^ and 

Tr (p ; ) entries, are constant in time. This should 

not be surprising as the effective Hamiltonian does not 
allow for the individual state populations to change. 

When just the two qubits are considered (trace over 
the cavity), the entanglement between them will undergo 
fluctuations. In line with the monogamy of entanglement 
every time when the cavity almost completely entangles 
to the qubits, the qubits themselves must share very little 
bipartite entanglement, because now the system forms a 
tripartite entangled state as a whole. 

When we calculate the amount of tripartite-shared en- 
tanglement we find 



Neg (Q1CIQ2) {t) = Neg (Q2CIQ1) (t) = ^ • (15) 

Note that the initial state of the cavity a has no effect 
on the results. It is important to note that under dis- 
sipationless evolution the entanglement between the two 
subsystems spanned by the joint qubit-cavity subspace 
and the other qubit does not change with time (i.e. there 
will be no bipartite entanglement variation between the 
two qubits), thus since Neg ((3i(52|C) > the total en- 
tanglement can only increase relative to its initial value. 

The behaviour of Neg (Q 1(52! C*) (see Fig. 2) displays 
periods of entanglement and disentanglement between 
the qubits and the cavity. This is due to the fact that ev- 
ery period of length 2n/5, the coherent state of the cavity 
returns to its initial state. Figure 2 also shows that the 
strength of qubits-cavity entanglement formed depends 
on detuning. The coherent states under detuned driving 
of the qubits change their amplitudes to a limited extent. 
The values of the coherent state amplitude and phase fol- 
low a circular trajectory in a complex plane centered at 
tto ± g/S with periods 5 and radii g/5, where ao is the 
initial coherent state amplitude. 

The creation of entanglement between the qubits and 
the cavity, however bares consequences to the qubit-qubit 
subsystem. Previously in Ref. [15, 16] we saw that qubits 
can undergo oscillations in their relative entanglement 
strengths (even if we take the k — ?> limit of equation 
(10)). By considering the solutions and Figure 2 we can 
see that throughout the evolution the qubit-qubit sub- 
system oscillates between completely entangled and par- 



C(c\Q^Q2) = W 1 - exp ^2 j , 

which incidentally is just twice the negativity in Eq. (14) 
[28]. If we calculate the C^cQi) {C(CQ2)) by taking the 
trace over the subsystem Q2 (Qi), we find them equal to 
zero. Thus we get that 

TCQ1Q2 = 1 - exp ( 1 

Thus we see that for small detunings most of the entan- 
glement is shared among the three entities and only when 
time t is close to an integer multiples of 2^/5 then the en- 
tanglement between the qubits subsystem and the cavity 
is lost, resulting in recovery of the bipartite entanglement 
between the qubits (see Figure 3). 

The total entanglement TE will simply be 
Neg ((5i(52|C') + 1, thus we can extend the conclu- 
sions above to the total entanglement in the system as 
their qualitative nature does not change. It is interesting 
to analyse the case when dissipation is present, which is 
the focus of the next section. 



V. DISSIPATIVE CAVITY 

Using the solutions (5), and repeating the analysis pre- 
sented above in the dissipative cavity case we find the 
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FIG. 3: (Color Online) A plot of a three-tangle and qubit-qubit subsystem concurrence for a dissipationless cavity (plum solid 
and green dashed lines respectively). We see that the at times of decreased qubit-qubit entanglement, the three-tangle grows, 
supporting our conjecture that the qubit-qubit entanglement is transferred to entanglement of a tripartite system. Plots made 
for 51 = 52 = 100 MHz, and S = 200 MHz (left) and S = 600MHz (right) 



negativities 

Neg(Qig2|C) = -1) 



frequency of re- and disentanglement. Additionally, we 
point out that of the qubit-qubit (non)-dissipative con- 
currences 



Neg(QiC|g2)(t) = Neg(g2C|gi)(t)--e''S 

with hi = hi (g, t) is given by Eg. (6) and A stems from 
the definition (7) and in this case reads 

2 _ / 4:g'^e^'^\cos{St) - cosh(Ki)) 

We can see that as a result of cavity dissipation the 
negativity Neg (giCig^) (t) (constant when k — be- 
comes a nontrivial function of time whose plots are pre- 
sented in Figure 4 for different values of dissipation rate 
K and detuning S. Other than the presence of detuning 
creating a (dis)entanglement oscillation, we see two com- 
peting effects playing a role here. Coupling-to-dissipation 
ratio at resonance leads to a decay of qubit-qubit en- 
tanglement and creation of qubits-cavity entanglement; 
detuning on the other hand, limits the qubit-qubit disen- 
tanglement, by means of impairing the qubit-cavity en- 
tanglement as we have seen in the previous section in 
Figure 2. Small detunings, facilitate formation of coher- 
ent states of larger maximum amplitudes which make the 
state resemble the GHZ state to a larger extend, thus dis- 
entangling the qubits (lowering the qubit-qubit concur- 
rence). When at the same time cavity dissipation rate is 
higher then the coherent state formed decreases its ampli- 
tude since a cx g/{K+i5) resulting in a slower monogamy 
of entanglement induced disentanglement rate. For large 
detunings, the qubits are coupled to the cavity field less, 
causing a smaller coherent state amplitudes and reduc- 
ing the disentanglement rate. Greater cavity decay rate 
again only amplify this process. This explains the find- 
ings of the previous paper [15, 16], and we see this pro- 
cess quite clearly here due to detuning which marks the 



C^{t) = C(0)exp 



4 (ffi + 92) 



{l- Kt 



it is the second one that decays slower in time. Thus 
the cavity dissipation process is not the cause for the 
qubit-qubit disentanglement, but rather it is the factor to 
decreases the initial disentanglement rate (Figure 4) . We 
can see this clearly, since the concurrence function C^, is 

always decreasing (^C^ — g^K~^^ with a long time limit 
equal to zero for any finite positive value of k, however for 
infinite k concurrence remains constant. This shows that 
in this system the large cavity decay prohibits the qubits 
from entangling to the cavity and protects the bipartite 
entanglement that the qubits share. 

To get a fuller picture of the negativity time evolution, 
we would need to find the three-tangle (where we use 
the K subscript to denote the dissipative case) which is 



1 
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The conclusion that we can draw from these results is 
that in either scenario there is entanglement being formed 
between the qubits and the cavity and as a consequence 
of monogamy of entanglement, the greater the degree of 
entanglement between the qubits and the cavity the less 
can they retain their inner-qubit entanglement. With the 
decay in the cavity present the entanglement between the 
qubits and the cavity weakens, which is a consequence of 
the cavity decay which can be seen as the cavity field 
entangling to the states of the environment which since 
traced over (a procedure carried out when deriving the 
Lindblad term) , result in a continuous degradation of any 
entanglement present in the system as a whole. 
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FIG. 4: (Color online) The plot of Neg{QiC\Q2) (t) (red-solid) and Neg(QiQ2|C) (t) (blue-dashed) made for gi + g2 ^ 200 
MHz, 5 = 100 MHz (first column), S = 300 MHz (second column) and 5 = 500 MHz (third column), and k = 100 MHz (first 
row) and k = 1 MHz (second row). Bottom-left: {g = 50 MHz, 5=1 GHz, k = 400 MHz) we see that the Neg {QiQ2\C) (t) is 
much longer lived and eventually decays but always later than Neg (Q1CIQ2) {t). 

We see that the entanglement generated between the cavity and the two qubits depends on a delicate balance between the 
detuning and the dissipation rate. Bottom-right: comparison between the qubit-qubit concurrence decay rate with and without 
decay. Plots made for gi + g2 = 200 MHz, n = {500, 100, 50, 0} MHz in blue, red, green, black (increasing dashing frequency) 
respectively. We see that larger cavity decay rates offer slower concurrence decays, however due to further cavity-environment 
coupling asymptotically concurrence is zero. 



VI. CONCLUSIONS 



In this article we have studied the dynamics of tri- 
partite entanglement between two driven qubits non- 
resonantly coupled to a cavity. Using tripartite entan- 
glement measures (negativity and three-tangle) we have 
shown that the previously reported entanglement loss 
followed by its revival is a consequence of an entangle- 
ment formation and subsequent disentanglement between 
the subsystem composed of qubits and the subsystem 
spanned by the cavity. Additionally, further tripartite 
entanglement loss is due to the dissipation of the cavity, 
which can be seen to form a greater entangled state with 
the environment states which has been traced over in a 
process of a derivation the master equation. From this 
one can see the qubits coupled to the cavity, which acts 



as an intermediate non-Markovian bath, which is further 
coupled to a Markovian one. This non-Markovian-like 
behaviour can be seen to be due to the presence of the 
external driving field and the cavity frequency mismatch 
S, which clocks the (dis)entanglement process. With this 
work we want to emphasize the danger of attributing 
all correlation losses to dissipation alone as seen by the 
evolution of correlations in this system, where the qubit- 
qubit entanglement is lost due to the qubits-cavity for- 
mation even for k = and it was only the formation of a 
larger system-environment entanglement formation that 
lead to tripartite intra-system entanglement loss. 
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FIG. 5: (Color online): Tangles/qubit-qubit concurrences, without (green dashed/plum solid) and with (blue dashed/red solid) 
dissipation present made for gi ~ g2 = 100 MHz, 5 — 100 MHz (first column), 5 = 500 MHz (second column), and k = 10 MHz 
(first row) and k, = 100 MHz (second row) We see that the entanglement generated between the cavity and the two qubits 
depends on a delicate balance between the detuning and the dissipation rate. If detuning is small the tangle is also very short 
lived. 
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[comment 1] We assume that such an entangled state is avail- 
able. In our follow up work we show how one could pro- 
duce an entangled state using a system of strongly driven 
qubits. 



